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ABSTRACT 

We present the formation of a Kinematically Decoupled Core (KDC) in an elliptical galaxy, resulting 
from a major merger simulation of two disk galaxies. We show that although the two progenitor 
galaxies are initially following a prograde orbit, strong reactive forces during the merger can cause a 
short-lived change of their orbital spin; the two progenitors follow a retrograde orbit right before their 
final coalescence. This results in a central kinematic decoupling and the formation of a large-scale 
(^2 kpc radius) counter-rotating core (CRC) at the center of the final elliptical-like merger remnant 
(M* = 1.3 x 10 11 Mq), while its outer parts keep the rotation direction of the initial orbital spin. 
The stellar velocity dispersion distribution of the merger remnant galaxy exhibits two symmetrical 
off-centered peaks, comparable to the observed “2-cr galaxies”. The KDC/CRC consists mainly of 
old, pre-merger population stars (older than 5 Gyr), remaining prominent in the center of the galaxy 
for more than 2 Gyr after the coalescence of its progenitors. Its properties are consistent with KDCs 
observed in massive elliptical galaxies. This new channel for the formation of KDCs from prograde 
mergers is in addition to previously known formation scenarios from retrograde mergers and can help 
towards explaining the substantial fraction of KDCs observed in early-type galaxies. 

Keywords: galaxies: elliptical and lenticular, cD — galaxies: formation — galaxies: interactions — 
galaxies: kinematics and dynamics — galaxies: stellar content — galaxies: structure 


1. INTRODUCTION 

Early-type galaxies (ETGs) are the end-products of 
complex assembly and evolutionary processes that de¬ 
termine their shape and dynamical structure. Signa¬ 
tures of such past processes in present-day ETGs are 
likely to be in the form of peculiar kinematic subsystems 
that reside in their central regions. Such subsystems are 
called Kinematically Decoupled Cores (KDCs) and they 
are defined as central stellar components with distinct 


galaxy (e.g. 

McDermid et al.12006 Krajnovic et al. 2011; 

Toloba et af 
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re first discovered using one-dimensional 
ctroscopic observations of the stellar kine- 

TGs dEfstathiou et al. 

1982 Bender 

1988; 

ingworth 1988). More recently, integral- 
ectroscopic surveys such as SAURON (Ba- 


con et al. 


CAL1EA ( Sanchez et al.|2012 ), being able to provide full 


(Cappellari et al. 2011), or 


two-dimensional observations 


of the stellar kinematics, 


have favored the detection of KDCs and revealed that 
a substantial fraction of ETGs in the nearby universe 
show kinematic decoupling in their central regions. This 
fraction ranges in different surveys, depending mainly on 
technical and sample-selection biases. 

Notably, the fraction of ETGs that host K DCs in 
the SAURON sample of 48 E+SO galaxies (de Zeeuw 
et al. [2002 ) is substantially high, especially in the cen- 
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slow-rotating ETGs: 8 out of the 12 slow rota- 


67%) fro m the main survey host a KDC (Em- 


sellem et al. 2007). In the ATLAS 3D volume-limited sam- 
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pie o f 260 ETGs, this fraction is 47% (Krajnovic et al. 
201 1|) . The KDCs found in slow rotators are typically 


“old and large”, with stellar populations that show little 
or no age differences with their h ost galaxy (older than 8 


Gyr) and sizes larger th an 1 kpc ( McDermid et al.| 2006 


Kuntschner et al. 2010). 


KDCs are also detected in fast-rotating ETGs: 25% of 
fast rotators from the main SAURON survey host KDCs. 
This type of KDCs are typically “young and compact”, 
with stellar populations younge r than 5 Gyr and size s 
less than a few hundred parsecs (McDermid et al.|[2006). 

We note that these fractions establish a lower limit 
to the true fraction of ETGs with kinematically decou¬ 
pled regions, considering projection effects, the fact that 
young and compact KD Cs are subject to techn ical or 
observational biases (e.g. McDermid et al. 2006), while 
many ETGs with resolved KDCs in their centers are sub¬ 
ject to different classificat ions throughout t he literature 
(e.g. 2cr-galaxies, see |Krajnovic et al.||2011). 

While a consensus is reached about the prominent exis¬ 
tence of KDCs in luminous ETGs, the physical processes 
and the rate at which they are formed are still poorly 
understood. Young and compact KDCs in fast rotators 
might have formed via star-formation in situ. Accord¬ 
ing to this scenario, the stellar component of the KDC is 
formed in initially kinematically misaligned gaseous re¬ 
gions, probably originating from externally accret ed gas 
or unequal mass merging (Hernquist & Barnes 1991), 
where the orientation of the merging orbit defines the 
orientation of rotatio n of the resulting KDC . Following 
this line of thought, Balcells & Quinn ( 1990|) suggested 
that counter-rotating cores can only result from retro¬ 
grade mergers. 

However, this scenario could not hold for the large and 
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old KDCs found in slow rotators, whose stellar popu¬ 
lation was probably formed at the same epoch as the 
main body of the galaxy. In this case, processes such as 
gas accretion or accretion of low-mass stellar systems are 
more likely to affect the outer parts of the galaxy and can 
not be consistent with observations that show no color 
gradi ents bet ween th e KDC and the surrounding galaxy 


1997). 


(Carolio et al. 

The most plausible formation scenario that could ex¬ 
plain the similarity of the stellar content of the KDC 
and the main body of the galaxy is major merging. 
This scenario has been confirmed in simulations (e.g. 


Bois et al. 2010, 2011), resulting in elliptical-like and 


slow-rotating merger remnants hosting KDCs only when 
the two progenitor galaxies were initially following ret¬ 
rograde merger orbits. However, observations indicate 
a lower limit to the true rate of occurrence of KDCs in 
ETGs which can not be explained only by retrograde 
mergers, pointing to the need of additional KDC forma¬ 
tion scenarios. 

Here we show that, a KDC might as well result from 
an initially prograde major merger. The kinematic de¬ 
coupling in the center of the final elliptical-like merger 
remnant can result from a short-lived change of the or¬ 
bital spin of the two progenitor galaxies right after their 
second encounter. This new channel for the formation of 
KDCs might serve as an additional mechanism that can 
help towards explaining their observed rate of occurrence 
in ETGs. 

2. SIMULATION PARAMETERS 


The simulation we use is described in IMoster et al.l 
(2011). It was perfo rmed using the TreeSPH-code 
GADGET-2 (|Springel||2005|), including star formation 
and supernova feedback. The two progenitor disk galax¬ 
ies are identical and they are composed of a cold gaseous 
disk, a stellar disk and a stellar bulge, which are embed¬ 
ded in a dark-matter and hot-gas halo. 

The gaseous and the stellar disk of each progenitor 
galaxy have exponential surface brightness profiles and 
they are rotationally supported, while the spherical stel¬ 
lar bulge follows a Hernqnist (1990) profile, and is ini¬ 
tially non-rotating The dark matter halo has a Hern- 


quist profile and a sp in parameter cons istent with cosmo 


logical simulations (Maccio et^al. 2008) . The hot gaseous 
halo follows the /3 -profile (Cavaliere & Eusco-Femiano 


1976|) and is rotating a round the spin axis of the disk 


unbound prograde parabolic orbit, with an eccentricity 
of e = 0.95 and a pericentric distance of r Pl = 13.6 kpc. 
Such an orbit is representative for t he most common ma- 


jor mergers in ACDM cosmology (Khochfar & Burkert 
2006). The two galaxies have an initial separation "of 
d s tart = 250 kpc. The orbital and the rotation spin of 
the first galaxy are aligned, while the spin axis of the 
second galaxy is inclined by 0=30° with respect to the 
orbital plane. The simulation lasts for 5 Gyr, such that 
the remnant elliptical galaxy is fully relaxed. 

3. MERGER REMNANT 
3.1. Structure of the merger remnant 

In order to connect the orbital and mass distribution 
of our simulated galaxy with observable properties, we 
create two-dimensional mock stellar mass maps as fol¬ 
lows. Stellar particles are projected such that the galaxy 
is seen edge-on with respect to the initial orbital plane of 
the merger. Particles are then binned on a regular grid 
centered on the baryonic center of mass of the galaxy. 
We adopt a distance of 20 Mpc, so that 1 arcsec corre¬ 
sponds roughly to 0.1 kpc. Our grid has a size of 20x20 
kpc, covering approximately twice the half-mass radius 
(r/i) of our galaxy, and a pixel size of 0.075 kpc, so that it 
corresponds to the spatial resolution of t he IRAC camera 


of the Spitzer Space Telescope (Fazio et al. 


We parametrize the galaxy’s projected ste. 


2004) 


lar mass dis- 


tribntion using the Multi-Gaussian Expansion (MGE) 


model (Monnet et al.|l992 


. 1992 | lEms ellem et alj!994|), as im¬ 
plemented by Cappellari~( 2002[ ). The intrinsic shape of 


the remnant’s stellar particle distribution is parametrized 
using an iterative method to obtain the best fitting ellip¬ 
soid to the distribution and to extract th e eig e nvalues o f 
the mass tensor inside this ellipsoid (Maccio et al. 2008). 
The intermediate-to-long and short-to-long axes ratios 
that we retrieve are p=0.88 and q=0.54, respectively, 
and the average projected ellipticity of the remnant is 
e = 0.49, estimated within 2 rv 

3.2. Kinematics of the merger remnant 

Stellar particles are projected along the chosen viewing 
angle and binned on a regular 20x20 kpc grid centered 
on the baryonic center of mass of the galaxy. In order 
to mimic real integral-field spectroscopic data, the pixel 
size of 0.1 kpc corresponds, at the adopted distance of 
20 Mpc, approximately to the spatial reso lution of the 


(see |Moster et al 2011, for a more detailed description SAURON spectrograph (Bacon et al. 2001). The bulk 


of the galaxy model). 

The stellar mass of each progenitor is M * = 5 x 10 10 M o 
and the bulge-to-disk ratio was chosen to be B/D = 
0.22. The mass of the cold gaseous disk is Mg, co id = 
1.2 x 10 10 M o , such that the gas fraction in the disk is 
23%. The virial mass of the dark matter halo is = 
1.1 x 1O 12 M 0 , while the mass of the hot gaseous halo 
is M g ,hot = IT x 1O 11 M 0 . The softening length is 100 
pc for stellar, 400 pc for dark matter and 140 pc for gas 
particles. 

The two progenitors are initially employed in a nearly 

1 We note that these initial structural properties of the progen¬ 
itors are influenced by their close interaction, i.e. they develop 
bulge rotation, bars and spiral arms in the first few hundred Myr 
of the simulation. 


velocity of the galaxy is estimated within a sphere of 
50 kpc around the center and subtracted from all parti¬ 
cle velocities. Then we extract the mass-weighted stel¬ 
lar line-of-sight mean velocity and velocity dispersion for 
every pixel. The extracted kinematic maps a re spatiall y 
binned using t he 2D Voronoi binning method (Cappellari 
fe Copin|20 03), based on a minimum number of particles 
per pixelTn the map. Signal corresponds to the number 
of particles per pixel and we adopt Poisson noise, such 
that our signal-to-noise ratio per bin ( SNbi n ) should cor¬ 
respond approximately to a target value SNt ~ 30. 

We also use a simple logarithmic fun ction inferred from 
CALIFA data (Husemann et al. 2013]) to construct mock 
velocity errors oFour binned kinematic data: 


Sv = hSN T (l + 1.4log N pix )/SNfon , kms 


-l 


(1) 
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Figure 1. From left to right: Component analysis of the merger remnant’s stellar population: all stars (old+young), old stars (>5 Gyr), 
young stars (<5 Gyr), and “Pr.Bulge” stellar population. The latter are particles that initially formed the bulge of the progenitor galaxies. 
Top row shows the 2D maps of the line-of-sight projected stellar mass surface density (S) and bottom row shows velocity dispersion (cr) 
for every component. Overplotted are the MGE best-fitting contours of the total projected stellar mass, as a reference. The remnant is 
seen edge-on with respect to the orbital plane of the two progenitors. 
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Figure 2. Line-of-sight projected stellar rotation for every component of Figure [T| The top row shows the 2D maps of the line-of-sight 
projected stellar mean velocity (v) for every component. The bottom row shows the stellar rotation curve extracted from the 2D velocity 
maps: black points correspond to the mean velocity per bin inside a slit of 1 kpc width along the apparent major photometric axis (the 
position of the slit is shown at the top row figures of Figure [TJ. Blue points correspond to bins outside the slit. Typical average errors are 
shown at the bottom. 
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p, a, p 2 a 2 p 3 p 4 



Figure 3. Separation (d) and specific orbital angular momentum 
(l z ) of the two progenitor galaxies as a function of time. Times 
Pi, P 2 , P 3 and p 4 indicate the time of the 1st, 2nd, 3d and 4th 
pericentric passage. Times a± and a .2 denote the 1st and the 2nd 
apocentric passage. The orbital spin of the two progenitors changes 
after p 2 . 


where N p i x is the number of pixels per bin. 

For the purpose of this work, we divide the stellar par¬ 
ticles of the remnant into 4 different components: “old 
stars”, which are stars that initially were part of the pro¬ 
genitors’ stellar material, (ages > 5 Gyr), “young stars”, 
which were formed during the merger (ages < 5 Gyr) 
and “all stars”, which is the total stellar content of the 
merger remnant. We also track the stellar particles in 
the remnant that initially formed the bulges of the two 
progenitor galaxies. These particles form the “Pr.Bulge” 
(Progenitor bulge) stars. The projected two-dimensional 
stellar mass and the stellar kinematics for every stellar 
component are shown in Figures [l] and [2| 

One can clearly see from the velocity maps the pres¬ 
ence of a large-scale KDC of radius ^2 kpc in the center 
of the elliptical merger remnant (Figure |2|. This compo¬ 
nent is a counter-rotating component with respect to the 
outer body of the galaxy and is most prominent in the 
“old” stellar population kinematics: stars that initially 
belonged to the two progenitor galaxies. 

On the other hand, “young” stars form a stellar disk 
which is almost aligned with the orbital plane of the two 
progenitors. This young stellar disk is rotationally sup¬ 
ported and strongly prograde-rotating, with a maximum 
velocity which is 4 times higher than the one of the “old 
stars”. We also note a weak sign of counter-rotation in 
the central region of the disk, seen in the extracted stellar 
rotation curve. 

Notably, stars that were initially part of the two 
progenitors’ bulges (“Pr.Bulge”) are globally counter¬ 
rotating in the merger remnant, exhibiting almost a 
solid-body rotation. 

One can also see the presence of two symmetrical off- 
centered peaks in the “all stars” stellar velocity disper¬ 


sion map of Figure [l] This feature is commonly observed 
in ETGs with counter-rotating com ponents (CRC) 


These galaxies are called “2cr-galaxies” (Krajnovic et al. 


2011), and they were ass ociated with exte rnal accretion 


of counter-rotati ng gas (|Rubin et a l. 1992) or major ret 
rograde mergers (Crocker et al. 2009f. Here we see that a 
2cr-galaxy results from a single, prograde major merger. 
We also note that the 2cr-feature is more prominent in 
the “all stars” map, where the young disk of stars is 
present, and less strong in the “old stars” map, even 
though the CRC is more prominent in the latter. In 
real ETGs a 2cr-feature usually implicates the existence 
of a CRC, while the opposite is more ambiguous. Here 
we show that the 2cr-feature might arise because of the 
presence of the CRC, but it is enhanced only if one of 
the components is fast-rotating 


4. ORIGIN OF THE KINEMATIC DECOUPLING 

In order to understand the origin of the kinematic de¬ 
coupling in the central region of the galaxy, we study 
the behavior of the merging orbits of its two progenitors: 
Figure [3] shows the separation (d) and the specific orbital 
angular momentum ( l z ) for one of the two progenitors as 
a function of time. The merging orbits are shown in Fig¬ 
ure |4j viewed face-on with respect to the initial orbital 
plane. 

At the time the two progenitors reach their first peri- 
center (pi=0.78 Gyr), they become tidally distorted, re¬ 
sulting in long trailing arms that expel loosely bound ma¬ 
terial from their disks. The orbital angular momentum is 
decreasing due to mass loss and dynamical friction (Fig¬ 
ure J3|. This causes the galaxies to approach their second 
pericenter (p2=2.08 Gyr) with almost radial orbits (Fig¬ 
ure [4|. After their second encounter, the two galaxies 
change their orbital spin and follow a retrograde orbit 
for ^300 Myr, until they finally reach their coalescence 
at t^2.4 Gyr. 

The sudden change from a prograde to a retrograde 
merger can be understood in the framework of reac¬ 
tive forces. Due to strong tidal interactions during the 
merger, the two progenitor galaxies are systems of vari¬ 
able mass; mass is constantly ejected along their short¬ 
lived trailing arms. The mass loss from each system 
results in a reactive force, known as the |Mestschersky| 
([1902} force : 

^ = m(V - 7) (2) 


where m is the mass loss rate, V the bulk velocity of the 
system and if is the velocity of the outflowing matter. 
The Mestschersky force acts upon the two galaxies as a 
“reactive thrust” which, if strong enough, can cause the 
change of the orbital spin. 

Figure [5] shows in detail this effect after the second 
pericentric passage. After the two progenitors approach 
closely, strong tidal forces that act upon them result in 
short-lived trailing arms, which mainly consist of their 
disks’ stellar component. Loosely bound material gets 
ejected along these arms, resulting into a strong reactive 
thrust to the main bodies of the progenitors, causing 
them to change their orbital spin and follow retrograde 


2 We should note, however, that most 2cr-galaxies do not exhibit 
a centrally peaked velocity dispersion, like the one presented here. 
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Figure 4. Merging orbits of the two progenitor galaxies, with the orbital plane seen face-on. Each curve (red and blue) corresponds to 
the orbit of each progenitor’s bulge center of mass; the red curve denotes the galaxy that has its disk aligned with the orbital plane. Red 
and blue arrows show the direction of the orbital spin. Left: pi denotes the 1st pericentric passage, when the two galaxies follow prograde 
orbits. The orbital spin changes sign right after the second pericentric passage P2 . Right: Zoom-in view of the region inside the dashed 
line on the right, when the orbits become retrograde. 


trajectories until their coalescence at t~2.4 GyfO 

The central region of the final remnant is counter¬ 
rotating and the width of the last oscillation before co¬ 
alescence corresponds to the size of the KDC ( rsj 2 kpc), 
which is prominent in the center of the galaxy for more 
than 2 Gyr after the kinematic decoupling of its progen¬ 
itors. 

Under this framework, one can explain why stars 
that were initially part of the progenitors’ bulges show 
global counter-rotation in the post-merger kinematics 
(Figure [2|: these stars, more tightly bound in the centers 
of the galaxies during their close encounters, can track 
the behavior of the orbital spin of their progenitors’ cen¬ 
ter of mass before coalescence. 

On the other hand, the outer parts of the galaxy keep 
the initial prograde spin. Gas and stars ejected during 
the merger are subsequently re-accreted, while inheriting 
the outer prograde spin, forming the prograde-rotating 
outer part of the remnant. 

We suggest that the Mestschersky force is present in 
every stage of the merger. We interpret the change of sign 
of the orbital angular momentum near the first apocenter 
ol i, as a result of this force (Figure [3]), which is also seen 
as a change of curvature of the two merging orbits near 
aq in Figure [4| However, at this time in the merging 
process the effect is not strong enough to change the 
orbital spin. 


5. SUMMARY AND DISCUSSION 

We have shown that a KDC in an early-type galaxy 
can result from an initially prograde major merger of 
two disk galaxies. This finding is in contrast to the com- 

3 The suggested mechanism could be responsible for the forma¬ 
tion of KDCs in non-retrograde close encounters, e.g|Baraes|(|2002|). 


monly suggested idea that KDC formation can only re¬ 
sult from retrograde mergers. We show the plausibility 
of an orbital reversal of a prograde merger, caused by 
reactive forces that act upon the two progenitors due to 
mass loss, which results in KDC formation in the final 
merger remnant. 

The KDC that resides in the center of the remnant 
shows strong counter-rotation for more than 2 Gyr af¬ 
ter the final coalescence of its progenitors. The KDC 
is most prominent in the old stellar population of the 
galaxy (ages > 5 Gyr) and is large in size (2 kpc radius), 
making it c onsistent with observat ions of KDCs in mas¬ 
sive ETGs (McDermid et al. 2006 ) and comparable t o 
the observed CRC/2-q galaxies ( Krajnovic et ah||2011 ). 
The fact that it results from an initially prograde merger 
provides a new channel for KDC formation that can add 
to the predicted rate of occurrence of KDCs and help 
towards explaining their observed fraction in ETGs. 

The suggested formation scenario depends on reactive 
(Mestschersky) forces, that act upon the progenitors due 
to mass loss during the merger, causing the reversal of the 
orbital spin. Since prograde mergers result in substantial 


mass loss compared to re trograde mergers (Toomre & 


Toomre 1972 Barnes 1988), we expect that such an effect 


is more likely to occur in prograde mergers. We would 
also expect this effect to depend on the mass ratio, the 
initial inclination, as well as the structural properties of 
the progenitor galaxies. Resolution effects might also 
influence properties o f the KDC^ such as its size and its 
position angle (e.g. |Bois et al.|[2010). 

Using the merger simulations presented 
(2011), we note that the formation of the 
depend on the particular form of feedback used or the 
specific values of hot and cold gas fractions employed in 
the progenitor galaxies. 


. in M oster et al. 
3 KDC does not 
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Figure 5. Stellar surface density of the two progenitor galaxies after their second pericentric passage (2.08 Gyr) where the orbital spin 
changes sign. During this time, strong tidal interactions acting upon the two progenitors result in short-lived trailing spiral arms (2.08-2.18 
Gyr), along which loosely bound material gets ejected from their main bodies. This results into a strong reactive force that causes the two 
galaxies to change their orbital spin and follow retrograde trajectories until their final coalescence at t~ 2.4 Gyr. The width of the last 
oscill ation before coalescence corresponds to the size of the KDC in the final merger remnant. Images were created with SPLASH fPrice] 
[2007] ). 


A larger statistic of merger simulations will allow us 
to understand how common this new channel for KDC 
formation is, which we plan to investigate in future work. 
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